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ABSTRACT 

The  strain  gauge  system  implemented  in  one  of  the  ROBOT-X  can¬ 
ards  is  described.  The  system  is  designed  to  measure  lift,  drag  and 
torque.  The  description  includes  design  considerations,  calibration 
procedures  and  data  reduction  techniques.  Data  are  presented  from  an 
early  ROBOT-X  flight  test. 
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DEFENCE  RESEARCH  ESTABLISHMENT  SUFFIELD 
RALSTON,  ALBERTA 


SUFFIELD  MEMORANDUM  NO.  1191 


IN-FLIGHT  LOAD  MEASUREMENTS  OF  THE  ROBOT-X  CANARDS  (U) 


by 


S.  Penzes 


1.0  INTRODUCTION 


ROBOT-X  is  a  rocket-boosted  aerial  target  designed  to  simulate 
a  variety  of  air  threats  for  use  in  exercising  air  defence  weapon 
systems.  Its  configuration  is  illustrated  in  Figures  1,  2  and  3.  The 
airframe  is  made  almost  entirely  of  composite  materials,  with  a  length 
of  3.4  m  and  wing  span  of  2.4  ra.  Propulsion  is  by  means  of  nineteen 
CRV-7  rocket  motors,  fired  in  pre-programmed  stages.  A  wide  variety  of 
flight  profiles  may  be  achieved  via  a  programmable  microprocessor-based 
digital  autopilot.  The  target  is  recovered  by  a  parachute  system. 
Flight  tests  conducted  during  1986  have  established  the  technical 
soundness  of  the  ROBOT-X  concept. 

The  ROBOT-X  airframe  is  unique  in  terms  of  both  aerodynamic 
configuration  and  structural  design.  Consequently,  there  is  consider¬ 
able  research  value  in  obtaining  full  scale  measurements  of  aerodynamic 
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loads  on  key  members.  The  canards  are  a  natural  choice  for  such  mea¬ 
surements,  because  they  are  of  primary  importance  from  a  control  stand¬ 
point  (pitch  and  altitude)  and  represent  one  of  the  most  unique  aspects 
of  the  ROBOT-X  configuration. 

This  memorandum  describes  the  strain  gauge  system  implemented 
in  one  of  the  ROBOT-X  canards  to  measure  lift,  drag  and  torque.  The 
description  includes  design  considerations,  calibration  procedures  and 
data  reduction  techniques.  Data  are  presented  from  an  early  ROBOT-X 
flight  test. 


2.0  DESIGN  APPROACH 


The  objective  is  to  produce  relevent  aerodynamic  data  for  the 
canards  from  data  acquired  during  the  f±  ghts  of  ROBOT-X.  This 
objective  was  accomplished  in  four  stages. 

1)  The  design  of  the  on  board  measurement  system.  This  design 
utilized  three  decoupled  strain  gauge  bridges.  These 
bridges  produce  output  which  is  directly  proportional  to 
the  instantaneous  lift,  drag  and  torque  on  the  canard. 

2)  The  determination  of  a  cross  talk  correlation  matrix  from 
preflight  testing  of  the  canard  with  known  loads.  This 
matrix  is  used  to  mathematically  remove  the  influence  of 
cross  channel  systematic  errors  from  the  inflight  data. 

3)  The  implementation  of  a  preflight  shunt  calibration  method¬ 
ology.  This  calibration  is  used  to  electronically  simulate 
loads  on  the  canard.  The  primary  purpose  is  to  verify  that 
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the  on  board  signal  conditioning  is  performing  as 
designed. 

4)  The  implementation  of  a  post  flight  analysis  method  that 
converts  the  flight  induced  strain  measurements  to  meaning¬ 
ful  aerodynamic  loading. 


3.0 


STRAIN  GAUGE  CONFIGURATION 


The  canard  configuration  and  internal  structure  are  shown  in 
Figure  4.  The  main  structural  member  is  an  aluminum  shaft  at  the  25% 
chord  line.  The  shaft  is  set  in  an  instant  set  polymer  (ISP)  casting, 
from  which  pockets  are  machined  to  reduce  the  weight  and  rotational 
inertia.  The  casting  is  covered  with  fibreglass. 

In  keeping  with  normal  aeronautical  practice,  the  principal 
loads  acting  on  the  canard  are  considered  to  be  lift,  drag  and  pitching 
moment.  Measurement  of  the  strains  produced  by  these  loads  is  accom¬ 
plished  by  mounting  strain  gauges  on  the  shaft  as  shown  in  Figures  5 
and  6.  This  particular  configuration  was  chosen  for  the  following  rea¬ 
sons. 


1)  If  gauges  // 1  and  //2,  used  to  measure  lift,  are  wired  into 
adjacent  arms  of  the  Wheatstone  bridge  (see  Figure  7) 
then: 

the  effects  of  axial  load,  as  would  be  introduced  by 
side  forces  due  to  high  roll  rates,  are  cancelled; 

-  the  effects  of  thermal  coefficient  mismatch  between  the 
shaft  and  the  strain  gauges  are  cancelled; 
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-  the  output  of  the  bridge  is  twice  that  of  a  single 
active  arm  bridge;  and 

-  the  effects  of  temperature  induced  strains  are  cancelled 
provided  it  is  assumed  that  the  temperature  at  the  cross 
section  is  uniform. 

2)  The  same  arguments  hold  for  the  drag  measurement  bridge 
(gauges  #3  and  it 4) . 

3)  The  torque  bridge  is  also  thermally  compensated  and  pro¬ 
vides  four  times  the  output  of  a  single  active  arm  bridge. 

With  the  gauges  mounted  on  the  shaft  as  shown  in  Figure  6  and 
considering  that  the  whole  assembly  (canard  surface  and  shaft)  rotates, 
it  can  be  said  that  this  configuration  provides  information  in  the 
"canard"  co-ordinate  system.  Conversion  to  the  "body"  coordinate 
system  is  easily  accomplished  through  a  coordinate  transformation  based 
on  the  sweepback  angle  of  the  canard  (8g)  ,  the  canard  incidence  angle 
(ac,  measured  with  a  potentiometer  in  the  servo  motor)  and  vehicle 
flight  path  angle  (/,  estimated  by  the  autopilot).  This  procedure  is 
described  in  Section  6. 


4.0  STATIC  CALIBRATION  PROCEDURE 

If  the  measured  quantities  of  strain  are  to  be  used  to  derive 
an  estimate  of  the  aerodynamic  forces,  consideration  must  be  given  to 
the  systematic  errors  due  to  factors  such  as: 
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gauge  installation  misalignment, 

the  fact  that  strain  gauges  actually  provide  an  average 
strain  over  the  active  grid  area,  and 

the  strain  gradient  produced  by  the  curvature  of  the  sub¬ 
strate  to  which  the  gauges  are  applied. 

It  was  determined  that  a  calibration  procedure  should  be 
established  which,  combined  with  some  mathematical  compensation,  would 
provide  a  means  to  compensate  for  these  systematic  errors. 

A. 1  Test  Set-Up  and  Procedure 

It  was  determined  that  the  best  method  to  perform  the  calibra¬ 
tion  was  to  place  the  canard  in  a  fixture  resembling  the  one  used  in 
ROBOT-X  and  then  to  load  the  canard  at  various  locations  while  monitor¬ 
ing  the  bridge  outputs.  This  canard  calibration  procedure  is  hereafter 
referred  to  as  the  CCP.  The  schematic  illustrating  the  test  set  up  is 
shown  in  Figure  8.  This  was  accomplished  by  using  a  milling  machine  as 
illustrated  in  Figures  9  to  11. 

The  loading  points  were  constrained  by  the  available  movement 
of  the  milling  machine  bed  while  the  magnitude  of  the  load  was  limited 
by  the  location  of  the  load  point.  This  last  limitation  was  based  on 
the  internal  geometry  of  the  canard.  As  noted  earlier,  pockets  are 
machined  out  of  the  canard's  ISP  casting  to  reduce  weight  and  rota¬ 
tional  inertia,  and  then  the  whole  casting  is  covered  with  fibreglass. 
If  a  load  point  falls  in  the  middle  of  an  air  pocket  (i.e.,  between 
ribs  and  between  the  25%  chord  line  and  the  trailing  edge)  the  entire 
load  would  be  supported  by  the  fibreglass  skin  alone.  Considering  that 
the  load  was  to  be  applied  over  a  0.25"  diameter  eraser  tipped  point 
(so  as  to  accommodate  skin  curvature)  a  practical  load  limit  of  20  lb 
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was  established.  This  translates  to  a  localized  pressure  of  approxi¬ 
mately  400  psi.  The  exception  to  this  limit  was  loads  applied  along 
the  25%  chord  line  where  the  skin  is  immediately  backed  by  ISP.  The 
allowable  load  limit  was  increased  to  30  lb  along  this  line  of  action. 
The  resultant  load  point  geometry  is  shown  in  Figure  12.  The  selection 
of  the  load  points  also  had  the  added  rationale  that  in  theory: 

a  load  at  point  #1  produces  an  output  on  the  bridge  mea¬ 
suring  drag  but  none  on  the  torque  or  lift  bridges, 
the  same  load  at  point  #4  produces  the  numerical  equivalent 
of  point  #1  but  on  the  bridge  measuring  lift.  There  will 
be  no  output  from  the  bridges  measuring  drag  or  torque, 
and 

the  same  load  applied  at  point  #7  produces  an  output  on  the 
bridge  measuring  torque,  the  same  lift  output  as  for  point 
// 4,  and  no  output  of  the  bridge  measuring  drag. 

4.2  Analysis  of  the  Canard  Calibration  Procedure 

4.2.1  Theoretical  Analysis 

a)  Lift  and  Drag 

o  My  My 

e  =  -  and  a  =  —  or  e  =  — 

E  I  El 

where  e  =  strain  (in/ in) 

o  =  normal  stress  (psi) 

E  =  Young's  Modulus  (10  x  10*  psi  for  aluminum) 

M  =  bending  moment  (in-lb) 


UNCLASSIFIED 


UNCLASSIFIED 


I  =  area  moment  of  inertia  (in*) 
y  =  distance  from  the  bending  neutral  axis  to  the 
point  for  which  the  stress  is  being  calculated. 


Substitution  of  the  appropriate  values  leads  to 


eL  =  eD  = 


where  F  =  applied  force 


d  =  moment  arm 


r  =  radius  of  the  shaft 


y  =  r  (maximum  stress  at  the  outermost  fibres  of  the 
shaft) . 


b)  Torque 


€„  =  —  and  x  =  —  or  e  =  - 

2G  J  2GJ 


where  x  =  shear  stress  (psi) 

G  =  shear  modulus  (psi) 

J  =  polar  moment  of  inertia  (in*) 

c  =  distance  from  the  shear  neutral  axis  to  the  point 
for  which  the  strain  is  being  calculated. 


However  from  [2] ,  a  circular  shaft  under  pure  torsion 
has  its  principal  strain  axes  rotated  45°  from  that  of  the 
shear  strain  axis.  The  gauges  pair  used  to  measure  the 
torque  on  the  shaft  are  tilted  at  ±45°  from  the  longitudi¬ 
nal  axis  (see  Figure  6)  and  thus  the  bridge  measures  the 
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principal  strains.  With  the  substitution  that  G 
equation  (2)  reduces  to 


E 

2(l+u) 


2Fi ( 1+u) 

eT -  (3) 

Enr3 


where 


c  =  r  (maximum  shear  at  the  outer  most  fibre  of  the 
shaft) 

8  =  torque  moment  arm 
F  =  applied  force 
u  =  Poisson's  ratio 


Strain  is  also  induced  by  a  shear  stress  effect  [3] .  These 
strains  are  however  cancelled  by  virtue  of  the  gauge  placements 
(Figure  5)  combined  with  the  placement  of  the  active  gauges  within  the 
wheatstone  bridge  (Figure  8).  For  example,  for  lift  [A]: 


AE 


out 


=  E. 
in 


(R1+RJ ) 3 


(A) 


where  E. 

in 


AE 


out 

AR, 


AR. 

*s 


=  voltage  applied  to  the  bridge 
=  voltage  output  of  the  bridge 
=  resistance  of  gauge  #1 

=  change  in  resistance  of  gauge  #1  due  to  effects 
other  than  shear  stress 

=  change  in  resistance  of  gauge  if  l  due  to  shear 
effects. 
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Mechanically,  R3  and  R3  are  mounted  on  opposite  sides  of  the 
shaft  and  therefore  produce  resistance  changes  of  opposite  signs.  For 
an  upward  force,  #1  is  in  compression  and  therefore  ARa  <  0  while  #2  is 
in  tension  and  ARa  >  0.  The  shear  effect,  on  the  other  hand,  is 
equivalent  to  the  load  divided  by  the  area  and  is  uniformly  distributed 
over  the  whole  section  (both  ARlg  and  AR2g  are  greater  than  zero). 

This  can  be  summarized  as: 

AR,  =  -AR,  and  AR,  =  AR, 

11  1s  *s 

With  the  added  provision  that  R3  =  R3  (which  is  reasonable  if  the  two 
gauges  are  from  the  same  lot  number) ,  substitution  of  these  values  into 
equation  (4)  yields 


E. 

in 


(5) 


Examination  of  equation  (5)  shows  the  independence  of  the  bridge 
outputs  from  the  shear  effect.  Equations  (1)  and  (3)  will  be  used  to 
determine  the  theoretical  strains  of  the  CCP.  Examination  of 
equations  (1)  and  (3)  shows  that  once  the  relevant  strain  is  known, 
each  equation  has  a  pair  of  unknowns.  These  unknowns  (F  and  d  or  i) 
can  be  grouped  as  an  unknown  moment  or  couple  so  as  to  make  these 
equations  determinate.  The  implementation  of  the  post  flight  analysis 
method  will  then  be  used  to  convert  the  three  moments  to  meaningful 
aerodynamic  loadings. 
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4.2.2  Experimental  Analysis 

The  analytically  determined  strains  must  undergo  an  additional 
correction  before  comparison  to  the  strains  measured  during  the  CCP  can 
take  place.  This  correction  must  be  done  to  compensate  for  two 
factors: 

1)  the  discrepency  between  the  actual  gauge  factor  and  the 
gauge  factor  set  on  the  portable  strain  indicator  (which 
was  set  to  give  a  direct  reading  of  the  applied  load) ; 

2)  the  correction  for  the  bridge  factor  introduced  when 
multiple  active  arm  Wheatstone  bridges  are  used. 

The  two  equations  of  concern  are: 

e.K.  =  e*K*  and  e!  *  e*(BF) 

XI  1  1  1 

analytically  determined  strain 
strain  corrected  for  gauge  factor 
desired  corrected  strain 

gauge  factor  on  the  portable  strain  indicator 
bridge  factor 

gauge  factor  for  strain  gauge  i 

then 


where 


e.  = 

l 

e*  = 
l 

e!  = 

l 

K*  = 
BF  = 
K.  = 

l 


e!  = 

l 


e.K.(BF) 

K* 


(6) 
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For  the  gauges  used  on  the  canards  for  flight  #1: 


=  2.15 


=  2.115 


K*  =  1.6 


BFl  =  BFd  =  2 


BFt  =  4 


Table  I  contains  a  summary  of  the  analytically  expected  strains 
(as  derived  from  equations  (1)  and  (3)  and  corrected  by  equation  (6)) 
as  a  function  of  load  magnitude  and  load  location.  Table  II  contains 
the  experimental  results  of  the  CCP. 


Comparison  of  Tables  I  and  II  shows  the  basic  accuracy  of  the 


theoretical  calculations.  The  minor  differences  can  be  attributed  to 


channel  crosstalk  induced  by  gauge  misalignments. 


4.2.3  Calculation  of  the  Crosstalk  Matrix.  [C] 


The  crosstalk  matrix,  [C] ,  is  required  so  as  to  provide  a  means 
to  convert  the  raw  measured  output  of  the  bridges  to  a  decoupled,  ana¬ 
lytically  relevant  set  of  data.  That  is 


eL  cn  c 1 2  c i 3  eL 

ep  ”  ^3  2  ^3  3  eD 

eT  ^31  ^32  ^33  eT 


or  [e1]  =  [C][e] 


The  data  obtained  from  the  CCP  is  ideally  suited  for  the  derivation  of 
[C]  since  [e]  is  the  measured  output  (Table  II)  and  [e1]  is  the  analy- 
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tically  corrected  strain  vector,  (Table  I).  Data  obtained  from  the 
most  extreme  locations  in  combination  with  the  highest  load  have  the 
highest  signal  to  noise  ratio.  Therefore,  load  points  3,  6  and  9  at  a 
load  of  20  lb  will  be  used  to  calculate  [C] .  It  will  be  shown  that 
only  one  triplet  is  required  to  derive  [C] . 


Load  point  3 


[e1]  = 


Load  point  6 


[€']  = 


0 

-5 

657 

and  [e]  = 

638 

0 

11 

657 

633 

0 

and  [e]  = 

A3 

0 

10 

_  __ 

_  _ 

Load  point  9 


[e*]  = 


— 

657 

639 

0 

and  [e]  = 

20 

-144 

-125 

_  _ 

The  determination  of  the  nine  elements  of  the  [C]  matrix  as 
indicated  in  equation  (7)  would  require  nine  independent  equations. 
This  is  possible  since  nine  load  points  were  used;  however,  if 
equation  (7)  is  inverted  so  as  to  read 


Ce]  =  [C]-i  [e'] 
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this  provides  a  simpler  means  to  determine  [C]  since  equation  (8)  can 
then  take  advantage  of  the  fact  that  the  majority  of  the  theoretically 


corrected  strains  are  zero. 


For  example,  if  equation  (8)  is  expanded  for  load  point  3 


cn  ci»  ci 


C21  C22  CJ1 


^31  ^33  ^33 


it  is  evident  that  the  elements  of  the  second  column  of  the  [C] * 1  mat¬ 
rix  can  be  uniquely  determined.  The  expansion  of  equation  (8)  for  load 
point  6  will  similarly  allow  the  determination  of  the  first  column  of 

[c]->. 


Because  the  applied  torque  is  coupled  to  a  generated  lift 
moment,  the  coefficients  defined  when  load  point  6  was  expanded  must  be 
utilized  in  conjunction  with  the  expansion  of  load  point  9,  to  derive 
the  third  column  of  the  [C] ~ 1  matrix.  The  result  is: 


0.963  -0.008  -0.044 


[C]*1  =  0.065  0.971  0.158 


0.015  0.017  0.936 
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[C] 


1.037  -0.067  -0.015 
0.008  1.032  -0.019 
0.047  -0.177  1.071 


In  an  ideally  decoupled  situation  the  [C]  matrix  would  be  an 
identity  matrix.  The  determined  [C]  matrix  compares  favorably  in  that 
its  major  diagonal  elements  are  near  to  unity  while  the  off-diagonal 
elements  are  small.  The  fact  that  the  off-diagonal  elements  are  not 
equal  to  zero  indicates  the  degree  of  coupling  between  any  given  input 
load  (lift,  drag,  or  torque)  and  the  expected  strain  measurements  (e^, 
eD  and  eT) . 

4.2.4  Performance  Checks  on  [C] 


The  worst  case  tested  is  at  the  most  inboard  load  points  with 
the  lowest  load.  As  a  test  case,  load  point  7  (at  10  lb)  will  be  ana¬ 
lyzed  and  the  results  will  be  compared  to  the  corrected  theoretical 
strains. 


[e]  = 


157 

5 

-65 
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from  which  [e‘]  =  [C]  [e]  = 


160 

6.2 

-72.1 


This  result  compares  favorably  with  the  results  in  Table  I  of 


[e*]  = 


164 

0 

-71.7 


Another  interesting  test  is  possible  if  the  principle  of  super¬ 
position  is  used.  Consider  the  following  scenario: 


10  lb  applied  at  load  point  3  while 
20  lb  applied  at  load  point  8  and 
30  lb  applied  at  load  point  4. 


This  results  in:  a  lift  moment  of  360  in-lb 
a  drag  moment  of  120  in-lb 
a  torque  of  -40  in-lb 


Ideally  the  measured  strains  would  be  the  sum  of  the  individual  strains 
from  Table  II,  that  is: 
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# 

I 


[e] 


1  +  479  +  473 

953 

315  +  12  +  33 

= 

360 

6  -  127  +  5 

-116 

_  _ 

_ 

Using  equation  (8)  with  the  determined  [C]"1  matrix  yields: 


[e] 

i 


986 

328 

-146 


Instead  of  comparing  these  results  to  the  strains,  they  will  be  con¬ 
verted  back  into  moments  using  the  reverse  process  of  equations  (6) , 
(3)  and  (1).  This  results  in  a  moment  vector  of: 


ml 

360.2 

md 

= 

119.8 

Mt 

-40.8 

_  _ 

_  _ 

This  compares  very  favorably  with  the  results  expected  from  the  princi¬ 
ple  of  superposition. 


The  determination  of  the  [C]  matrix  was  based  on  the  expansion 
of  loads  for  one  triplet  of  load  points.  This  in  turn  requires  that 
the  following  two  questions  be  answered: 
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1)  Does  [C]  vary  as  a  function  of  applied  load? 

2)  Does  [C]  vary  as  a  function  of  location? 


These  two  questions  can  best  be  answered  by  examining  the 
effects  of  applying  10,  20  and  30  lb  at  each  of  load  points  4,  5  and  6. 
The  comparison  will  be  done  by  converting  the  measured  strains  to  their 
respective  moments  and  then  comparing  the  moment  to  the  actual  moment 
that  would  be  generated  by  the  given  load.  This  comparison  will  be 
done  by  expressing  the  calculated  moment  as  a  percentage  of  the 
expected  moment.  Table  III  contains  the  resulting  data. 


Following  are  some  observations  of  Table  III; 


1)  The  normalized  lift  varies  from  98.0%  to  101.1%.  This 
variation  does  not  appear  to  have  a  trend  when  either  the 
load  points  are  kept  the  same  and  the  load  is  increased  or 
when  the  load  is  kept  the  same  and  the  application  point 


moves . 


2)  The  normalized  drag  varies  from  -0.6%  to  2.4%.  Again  there 
appears  to  be  no  trend  either  as  a  function  of  span  or  of 
load. 


3)  The  normalized  torque  varies  from  -0.4%  to  0.1%.  Again 
there  appears  to  be  no  trend  either  as  a  function  of  span 
or  of  load. 


4)  The  deviation  of  point  6  with  a  20  lb  load  should  be  100% 
lift,  0%  drag  and  0%  torque  since  this  point  was  used  in 
generating  the  [C]  matrix.  The  errors  can  be  attributed  to 
roundoff  errors  in  the  calculations. 
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The  conclusion  to  be  reached  from  the  examination  of  Table  III 
is  that  since  the  deviation  of  the  non-loaded  axis  (in  this  case  drag 
and  torque)  varied  only  from  -0.4%  to  2.4%  maximum  and  this  is  well 
within  engineering  limits,  if  any  spanwise  or  load  magnitude  dependency 
exists  in  the  canard  strain  measurements,  it  can  be  ignored. 


5.0  PREFLIGHT  CALIBRATION 


For  each  flight  of  R0B0T-X  it  was  decided  that  two  preflight 
checks  should  be  considered.  If  a  new  set  of  instrumented  canards  was 
to  be  installed  on  the  airframe  then  a  canard  calibration  procedure  was 
to  take  place.  The  intent  of  the  procedure  was  to  determine  the  proper 
crosstalk  coefficient  matrix  for  the  canard  in  question.  This  proce¬ 
dure  is  previously  outlined  in  Section  4.0. 

The  second  pre-flight  test  was  to  be  a  shunt  calibration  of  the 
strain  gauge  bridges.  This  procedure  would  be  done  on  the  assembled 
airframe  and  would  utilize  the  telemetry  subset  to  be  flown  with  the 
particular  canard.  This  procedure  would  provide  three  pieces  of  useful 
information.  First  of  all,  it  shows  that  the  outputs  of  the  bridges 
under  a  simulated  load  are  as  expected.  Secondly,  it  demonstrates  that 
the  telemetry  (signal  conditioning  plus  transmitter/receiver  plus  the 
real  time  processing  functions)  is  operating  properly.  Thirdly,  the 
calibration  points  that  are  required  by  the  telemetry  post-processor  in 
the  analysis  of  the  flight  test  data  are  determined. 

5.1  Shunt  Calibration 


Figure  13  shows  a  schematic  of  the  in-flight  measurement  system 
for  the  canard  strain  gauge  bridges.  As  a  terminal  strip  was  required 
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to  make  the  connection  between  the  gauges  on  the  canard  and  the  tele¬ 
metry  system,  it  was  decided  that  this  would  be  the  most  convenient 
location  to  introduce  the  shunt  calibration  resistors.  This  did,  how¬ 
ever,  require  that  the  lead  lengths  between  the  terminal  strip  and  the 
gauges  themselves  be  considered  in  the  calculation  of  the  shunt  resis¬ 
tor  value.  This  is  necessary  since  the  resistance  in  these  leads 
causes  an  overall  desensitization  of  the  bridge.  The  schematic  of 
shunt  calibration  is  shown  in  Figure  14.  The  relationship  between  the 
required  change  in  resistance  of  the  gauge  due  to  the  introduction  of 
the  shunt  resistor  and  the  equivalent  strain  is  [5] : 

AR 

KR 

g 

where  e  =  strain  to  be  simulated 

R  =  resistance  of  the  strain  gauge 

g 

AR  =  change  in  resistance  required  to  simulate  the  given  strain 
K  =  gauge  factor. 

However  K  must  be  rewritten  to  include  the  desensitization  due  to  lead 
resistance: 

R 

g 

K'  =  K  - 


where  R^  =  lead  resistance.  Therefore, 


UNCLASSIFIED 


AR 


Ve 


e  = 


or  e  * 


20. 


(10) 


K'R 


K'R 


where  R„  is  the  equivalent  resistance  of  the  parallel  circuit  of  R  and 
“  S 

Rg.  Equation  (10)  can  therefore  be  rewritten  as 


R  R 
s  g 


R  - 


R  +R 
s  g 


(ID 


K'R 


substituting  in  the  relationship  for  K'  and  solving  for  R  ,  yields 


R  (1-Ke)  +  Rl 


R  = 
s 


(12) 


Ke 


Equation  (12)  is  used  for  determining  the  shunt  calibration 
resistor  for  a  single  active  arm  wheatstone  bridge;  however,  the 
bridges  used  on  the  canard  are  multiple  active  arm  and  therefore  com¬ 
pensation  for  the  bridge  factor  must  be  introduced.  The  result  is: 


R  (1-Ke)  +  Rl 


R  = 
s 


(13) 


Ke  (BF) 


where  BF  =  bridge  factor  (1  for  single  active  gauge,  2  for  two  active 
gauges,  4  for  four  active  gauges). 

Since  a  shunt  calibration  can  only  produce  a  decrease  in  equiv¬ 
alent  resistance,  consideration  must  be  given  to  the  relationship 
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between  the  shunt  calibration  and  the  mechanical  aspects  of  the  strain 
gauge  bridges. 

The  following  force  sign  convention  will  be  used: 

1)  Positive  lift  causes  a  pitch-up  of  the  airframe. 

2)  Positive  drag  impedes  the  forward  velocity  of  the  air¬ 
frame. 

3)  Positive  torque  causes  pitch-up  of  the  airframe. 

With  this  convention,  positive  lift  would  be  out  of  the  page  in 
Figure  5  and  would  therefore  cause  gauge  //I  to  be  in  compression 
(AR  <  0)  and  gauge  #2  to  be  in  tension  (AR  >  0) .  This  therefore  leads 
to  the  conclusion  that  gauge  //I  is  to  be  shunt  calibrated  for  positive 
lift  and  gauge  #2  is  to  be  shunt  calibrated  for  negative  lift.  The 
same  logic  is  applied  to  the  drag  and  torque  bridges. 

The  determination  of  the  strain  value  to  be  used  in 
equation  (13)  is  based  upon  a  combination  of  material  sciences  plus  the 
flexibility  of  the  signal  conditioning  units. 

The  canard  shaft  is  manufactured  from  7075-T6  Aluminum  whose 
rated  strain  failure  is  7500  microstrain  and  therefore  the  max/min 
limits  were  set  to  ±7500  microstrain  for  the  first  flight.  Although  it 
was  estimated  that  the  canards  were  in  actuality  overdesigned,  it  was 
decided  that  a  complete  record  would  be  beneficial  if  failure  occurred. 
It  was  estimated  that  these  limits  would  be  refined  on  a  flight  to 
flight  basis. 

Having  requested  a  max/min  value  for  the  strain  bridges,  the 
designers  of  the  signal  conditioning  modules  then  attempted  to  provide 
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this  value  by  the  selection  of  appropriate  gain  control  resistors  with¬ 
in  the  signal  conditioning  modules . 

It  was  decided  that  the  most  convenient  location  to  do  the 
shunt  calibration  was  at  the  Building  15  Test  Complex  at  DUES  where  the 
ROBOT-X  vehicles  are  assembled  and  systems-integrated  and  which  is  also 
the  location  of  the  Data  Acquisition  Center  (DAC) . 

The  shunt  calibration  procedure  is  as  follows: 

1)  Assemble  all  the  components  to  be  used  for  a  given  flight. 
This  includes  the  appropriate  canard,  signal  conditioning 
units  and  telemetry  subset. 

2)  Place  the  avionics  beam  where  easy  access  to  the  signal 
conditioning  is  available  and  make  the  appropriate  connec¬ 
tions  to  the  canard  strain  gauges.  In  addition  a  hardwire 
connection  from  the  telemetry  transmitter  to  the  telemetry 
receiver  (in  the  DAC)  is  made. 

3)  Since  the  teleraetred  data  post-processing  software  (IDEAS) 
allows  the  use  of  multiple  calibration  points,  up  to  seven 
shunt  calibration  values  are  determined  for  each  bridge. 
These  calibration  points  are  used  to  convert  the  raw  data 
to  engineering  units. 

4)  The  telemetry  subset  is  powered  up  and  initial  bridge 
balances  are  done  for  the  strain  gauge  bridges  (zero  strain 
value) .  Since  the  telemetry  data  words  are  10  bits  wide, 
each  bridge  is  balanced  for  a  digital  count  of  512  (one- 
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half  of  210).  This  balance  is  observed  in  real-time  in  the 
DAC. 

5)  The  shunt  resistors  are  then  applied  and  the  appropriate 
change  in  the  bit  count  is  recorded  for  incorporation  into 
the  IDEAS  program  for  post  flight  analysis. 

6)  This  procedure  is  repeated  for  all  the  strain  gauge 
bridges. 

Incorporation  of  the  two  calibration  procedures  is  meant  to 
insure  that  not  only  is  the  in-flight  instrumentation  functioning  as 
designed  but  also  that  the  data  received  will  be  usable  in  the  analysis 
of  the  performance  of  the  vehicle. 

5 . 2  Structural  Re s onance 


After  the  first  flight  of  ROBOT-X  in  mid-March,  1986,  an  inter¬ 
esting  effect  was  noted  in  the  response  of  the  canard  strain  gauge 
bridges.  There  appeared  to  exist  a  resonance  in  the  output  of  the  drag 
bridge.  A  simple  measurement  of  the  wave  length  of  the  major  peaks 
showed  a  frequency  of  approximately  40  Hz.  In  an  effort  to  ascertain 
the  cause  of  the  resonance,  it  was  decided  that  a  rudimentary  struc¬ 
tural  analysis  would  take  place  on  the  next  airframe.  This  was  accom¬ 
plished  by  using  a  Hewlett  Packard  spectrum  analyzer.  Since  only  the 
frequency  was  of  major  concern,  the  bridges  were  hooked  directly  to  the 
analyzer  and  the  canard  was  excited  by  striking  it  a  sharp  blow  with  a 
rubber  mallet.  After  performing  a  Fast  Fourier  Analysis  the  results 
were  as  follows: 
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1)  The  lift  bridge  has  a  natural  frequency  of  approximately 
37  Hz. 

2)  The  drag  bridge  has  a  natural  frequency  of  approximately 
36  Hz. 

3)  The  torque  bridge  has  a  natural  frequency  of  approximately 
60  Hz. 

This  resonance  appeared  to  result  from  a  combination  of  the  structural 
response  of  the  canard  with  the  compliance  of  the  attachment  mechanism. 
It  was  also  noted  that  any  preload  in  the  axis  concerned  tended  to 
severely  limit  the  amplitude  of  the  oscillation  in  that  axis  while 
slightly  increasing  the  natural  frequency  and  only  moderately  decreas¬ 
ing  the  amplitude  in  the  orthogonal  axis.  Typically  the  resonant  fre¬ 
quency  of  the  orthogonal  axis  shifted  upwards  approximately  10  Hz. 

Fortunately,  a  simple  solution  was  available  for  removing  this 
undesirable  measurement  characteristic  from  the  data.  The  original 
purpose  for  installing  the  bridges  was  to  measure  the  strains  produced 
by  aerodynamic  forces  and  the  secondary  purpose  was  to  attempt  to 
ascertain  the  frequency  response  of  the  canards.  To  this  end,  the  low 
pass  filters  in  the  signal  conditioning  modules  were  tuned  to  100  Hz 
(the  estimated  frequency  response)  and  as  such  were  not  filtering  the 
40  Hz  structural  resonance.  The  answer  was  simply  to  retune  the  low 
pass  filters  to  10  Hz.  With  the  18  dB/octave  rolloff  of  the  filters, 
this  meant  that  the  structural  resonance  would  be  about  35  dB  down  at 
40  Hz.  The  validity  of  this  modification  was  proven  in  subsequent 
flights  when  no  undesirable  resonances  were  observed  in  the  strain 
gauge  data. 
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6 . 0  REDUCTION  OF  THE  IN-FLIGHT  DATA 


Having  taken  reasonable  precautions  to  ensure  that  the  data 
received  from  the  canards  are  valid,  the  method  to  analyze  the  data  had 
to  be  established.  The  output  of  the  strain  gauges  can  be  reduced  to 
only  a  moment  about  an  appropriate  axis  (see  equations  (1)  and  (3))  and 
therefore  whatever  the  shape  of  the  canard  loading  profile,  the  moments 
could  be  resolved  as  a  triad  of  forces  acting  at  some  point  on  the 
canard.  It  is  the  intent  of  the  analysis  to  solve  for  this  triad  of 
forces  and  its  application  point  on  the  canard. 

Figure  15  shows  the  relationship  of  the  force  triad  to  the 
global  coordinate  system.  Note  that  to  convert  the  strain  gauge 
coordinate  system  to  the  global  coordinate  system,  two  rotations  are 
required.  The  first  is  a  rotation  about  the  "z"  axis  through  an  angle 
of  0g  which  is  simply  the  sweepback  angle  of  the  X  chord  line  of  the 
canard  relative  to  the  fuselage  reference  line.  The  numerical  value  is 
derived  from  the  geometry  of  the  canard  configuration  and  is  39° .  The 
second  is  a  rotation  about  the  y'  axis  through  an  angle  of  0^  (which  is 
the  sum  of  y  +  ac) .  In  an  idealized  situation  the  unknown  force 
triplet  acts  at  an  application  point  (AP)  which  is  defined  as  the 
intercept  of  the  X  chord  line  of  the  canard  with  a  line  drawn 
perpendicular  to  the  X  chord  line  which  passes  through  the  area 
centroid  of  the  canard  planform  (see  Figure  16).  If  the  distance  from 
the  centre  of  the  strain  gauge  group  to  the  AP  is  called  "A",  then  in 
the  ideal  case,  the  coordinates  of  AP  in  the  x'y'z'  system  are: 


AP  =  (0,  -A,  0) 
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Since  the  canards  are  non-ideal  the  above  equation  requires  that  devia¬ 
tions  be  introduced  in  defining  the  location  of  the  AP.  Therefore,  let 
the  coordinates  of  AP  be  given  by: 

AP  =  (£,,  -A  -  is) 

where  £1#  and  £,  are  deviations  due  to  manufacturing  and  aerodynamic 
effects. 


Based  upon  the  BOC  drawings  and  physical  locations  of  the 
strain  gauges. 


A  =  9.928”. 


If  Rj ,  Rj  and  R3  are  the  magnitudes  of  the  force  triplet  in  the  xyz 
coordinate  system,  then  it  can  be  shown  that  a  rotation  about  "z" 
through  "0g"  and  then  a  rotation  about  "y"'  through  "0^"  yields 


R. ' 

COS0  COS0 

cos0  sin8 

-sin8 

1 

P  s 

P  s 

P 

V 

= 

-sin0 

s 

COS0 

s 

0 

V 

sin8  cos6 

P  s 

sin8  sin0 

P  s 

cos0 

P 

Ri 

R, 


(14) 


where  Rl',  R, '  and  Rs '  are  the  magnitudes  of  the  force  triplet  in  the 
x’y'z’  coordinate  system. 

Equation  (14)  is  the  desired  result  since  it  provides  the  tri- 
axial  components  of  the  unknown  force  triad  in  the  strain  gauge  coor¬ 
dinate  system. 
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yields 


Measuring  the  distances  to  the  CP  in  the  x'y'z' 


coordinates 


. — 

dx' 

V 

= 

-(A+£j) 

(15) 


Note  that  the  application  point  AP  is  assumed  to  rotate  about 
the  y'  axis  while  the  force  triplet  is  assumed  to  remain  in  the  global 
coordinate  system. 

The  strain  gauge  bridges  are  configured  to  measure  the  resul¬ 
tant  moments  and  therefore  the  moments  generated  by  [F]  must  be  found. 
This  process  is  the  simple  summation  of  moments.  Using  the  sign  con¬ 
vention  of  Section  5.1, 


M  . 
x' 

“  V 

+  <dZ-> 

Ra ' 

(16) 

V 

=  -Mx.)  V 

+  «*•> 

Ri' 

(17) 

“  (dx’}  V 

-  V 

V 

(18) 

The  force  components  required  to  determine  equations  (16),  (17) 
and  (18)  are  found  in  equation  ( 1 A)  while  the  displacement  components 
are  found  in  equation  (15).  An  examination  of  these  equations  shows 
that  there  will  be  six  unknowns:  RJt  R2 ,  RJf  £lt  and  however, 
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there  are  only  three  equations.  To  make  this  a  determinate  system  the 
following  assumptions  are  made: 


1)  R, 


0  This  states  that  the  side  force  F  is  zero  which  is 

y 

valid  under  typical  flight  conditions,  i.e.,  angle 
of  sideslip  is  approximately  zero. 


2)  E2  =  0  This  states  that  the  load  point  AP  does  not  vary 
along  the  "y"'  axis  which  is  also  valid  under 
typical  flight  conditions. 


3)  E,  =  0  This  states  that  the  AP  deviation  in  the  "z"'  axis 
is  zero  which  is  a  reasonable  assumption  due  to  the 
easily  manufactured,  thin,  symmetric  cross-section 
of  the  canard. 


Accepting  these  assumptions,  substituting  the  required  values 
from  equations  (14)  and  (15)  into  equations  (16),  (17)  and  (18)  and 
then  reducing  yields 


~ ~ 

M  , 

0  R.sin8  +R,cos0  -R.sinS 

x’ 

1  p  3  p  1  s 

M  , 

= 

-R,sin8  cos8  -R,cos9  0  R,cos8  cos©  -R,sin0 

y 

i  p  s  3  p  *ps*p 

M  , 

-R.sin9  R,cos8  cos0  -R,sin8  0 

1  s  1  p  s  3  p 

z' 

l 


i 


A 

0 


(19) 


Making  the  substitutions: 


K.  =  sin9  cos9 
1  p  s 

K,  =  cos8 
J  P 

K,  =  sin8 
5  s 
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K.  =  cos8  cos0 
4  p  s 

K,  =  sin0 
s  P 


into  equation  (19)  and  rewriting  yields: 


M  , 
x 

M  , 

y 

M  , 
z' 


A(R1K1  +  R3K2) 

-SjCRjK,  +  R3K2) 

~e1K3  +  A(R1Klt  -  R3Ks) 


Solving  these  three  equations  simultaneously  yields: 


e,  = 


-M  ,  A 

y 


m  , 

x' 


amv,k3  K5M  . 


M  ,  - 
z 


M  , 
x' 


K, 


R,  = 


K5K3 

A( -  +  KJ 

K, 


(20) 


(21) 


R,  = 


M 

x'  RjKj 


K2A  k2 


(22) 


Note  that  equation  (22)  requires  that  equation  (21)  be  solved  first. 


As  an  illustration  of  the  procedures  involved,  a  sample  calcu¬ 
lation  will  be  performed. 
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7.0  A  POST  FLIGHT  ANALYSIS  NUMERICAL  EXAMPLE 

As  an  illustration  of  the  procedure  to  be  used  in  the  reduction 
of  the  flight  test  data,  Figures  17  through  21  show  the  data  traces 
from  the  first  few  seconds  of  flight  //2  which  took  place  in  late 
July  1986.  The  launch  took  place  at  approximately  52.2  seconds  on  the 
time  scale  in  the  Figures.  These  graphs  were  produced  by  the  IDEAS 
post-processor  and  the  appropriate  calibration  constants  were  used  to 
convert  the  ordinate  of  each  graph  into  engineering  units.  As  an  arbi¬ 
trary  example,  consider  the  time  T  =  54  sec.  At  this  point: 

=  1575  pe  eD  =  48  pe  e^  =  -16  pe 
r  =  5.2°  ac  =  6.5° 

As  this  was  the  second  instrumented  canard,  a  pre-flight  CCP  was  per¬ 
formed  and  the  resultant  crosstalk  matrix  determined. 


[e] 


1.014  0.063  -0.053  1575 

-0.080  1.035  -0.054  48 


0.003  -0.154 


from  which  [e] 
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Using  equations  (1)  and  (3)  to  derive  the  moments  yields 

M.  =  1564  =  M  , 

L  x’ 

M  =  149  =  M  , 

D  z' 

Mt  =  -150.9  =  My, 

Substituting  =  39°,  0^  =  a ^  +  f  =  6.5  +  5.2  =  11.7°,  and  the  above 
derived  moments  into  equations  (19)  through  (22)  yields 

Sj  =  0.958" 

Rt  =  60.1  lb 
R3  =  151.2  lb 

This  result  is  illustrated  in  Figure  22.  At  that  instant  in  time 
(T  =  54  secs)  with  the  given  angle  of  attack  of  the  canard,  these 
derived  values  will  produce  the  strains  as  measured  by  the  strain  gauge 
bridges. 


The  purpose  of  this  example  is  to  simply  illustrate  the  proce¬ 
dure  used  to  derive  useful  aerodynamic  data.  A  detailed  analysis  to 
ascertain  the  significance  of  these  numbers  will  be  the  subject  of  a 
subsequent  report  which  will  be  issued  when  the  analysis  of  data  from  a 
number  of  flights  has  been  completed. 


8.0  CONCLUDING  REMARKS 


The  intent  of  this  report  has  been  to  illustrate  the  means  to 
measure  the  aerodynamic  forces  on  the  ROBOT-X  canards.  This  was  accom¬ 
plished  by  designing  a  transducer  circuit  composed  of  three  strain 
gauge  bridges.  As  for  any  transducer,  a  means  to  calibrate  and  verify 
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the  output  of  these  strain  gauge  bridges  was  developed.  To  date,  the 
general  trends  of  the  strain  data  correlate  well  with  what  can  be 
inferred  from  the  other  on-board  transducers,  the  autopilot  responses, 
and  the  performance  of  the  vehicle. 
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Lift 

Drag 

Torque 

Lift 

Drag 

Torque 

G  * 

L 

£'d 

e'T 

£'l 

£,d 

£'t 

0 
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0 

0 
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0 

0 
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0 
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0 
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0 
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0 
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0 
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0 

-144 

246 

0 

-71.7 

493 

0 

-144 

328 

0 

-71.7 
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Load  =  30  lb 


G  1  G  1  G  1 

L  D  T 


Units  are  microstrain. 
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Load  Point 


Load  =  10  lb 


Load  =  20  lb 


Load  =  30  lb 


iraber 

Fig.  12) 

Lift 
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Drag 
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Torque 

€t 

Lift 
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Drag 
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Torque 

eT 

Lift 

*L 

Drag 
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-3 
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4 
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33 

7 
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6 

317 

21 

5 

633 

43 

10 
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57 

7 
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5 
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10 
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- 
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8 
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6 
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12 
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9 
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8 
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20 
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TABLE  II  EXPERIMENTAL  RESULTS  OF  THE  CCP 
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Load  =  10  lb 

Load  =  20  lb 

Load  =  30  lb 

1 

Lift 

Drag 

Torque 

B 

Lift 

Drag 

Torque 

B 

Lift 

Drag 

Torque 

M 

■ 

58.8 

0.6 

-0.1 

120.6 

2.9 

-0.3 

179.4 

0.5 

-0.7 

1 

(98.0)* 

(1.0) 

(-0.2) 

60 

(100.5) 

(2.4) 

(-0.3) 

120 

(99.7) 

(0.3) 

(-0.4) 

180 

89.5 

0.01 

0.1 

179.4 

0.4 

-0.1 

270.4 

-1.1 

-0.3 

5 

(99.4) 

(0.01) 

(0.1) 

90 

(99.7) 

(0.2) 

(-0.1) 
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(100.1) 

(-0.4) 

(-0.1) 
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120.3 

-0.1 

0.03 

240.1 

0.1 

0.03 

364.0 

-2.3 

-2.9 

6 

(100.3) 

(-0.1) 

(0.03) 

120 

(100.0) 

(0.04) 

(0.01) 

240 

(101.1) 

(-0.6) 

(-0.8) 

360 

is  the  actual  value  of  the  lift  moment  generated,  which  is  used  in  normalizing 
the  measured  moments. 

*  Normalized  value  calculated  by  1/L^  times  measured  moment  (percent) . 
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Figure  3 

ROBOT -X  CONFIGURATION 
(Aft  View  Looking  Forward) 
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SECTION  A  -  A 


Figure  5 

LOCATION  OF  STRAIN  GAUGES 
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GRID  AREA  AT  45°  TO 
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TO  CANARD 


Figure  6 

GAUGE  ORIENTATIONS  ON  THE  CANARD  SHAFT 
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LIFT  MEASUREMENT 


DRAG  MEASUREMENT 


TORQUE  MEASUREMENT 


Figure  7 


NOTE: 


Rc  -  COMPLETION  RESISTOR 
R.  -  STRAIN  GAUGE  NO.  * 


BRIDGE  FACTOR  =  2 


BRIDGE  FACTOR  =  2 


7  BRIDGE  FACTOR  =  4 


WHEATSTONE  BRIDGE  CIRCUITS 
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Figure  9 

CCP  TEST  SET  UP 
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Figure  10 


OVERVIEW  OF  THE  CCP  TEST  INSTRUMENTATION 


LOAD  CELL  USED  DURING  THE  CCP 
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SCHEMATIC  OF  THE  IIM  -  FLIGHT  MEASUREMENT  SYSTEM 
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Figure  16 

THE  DETERMINATION  OF  THE  CONSTANT  "A" 
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Figure  17 

TRACE  OF  LIFT  FOR  FLIGHT  »2 
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